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We develop a theory for non-planar interaction between two identical type I spatial solitons propa- 
gating at opposite, but arbitrary transverse angles in quadratic nonlinear (or so-called x' 2 ' ) bulk media. 
We predict quantitatively the outwards escape angle, below which the solitons turn around and collide, 
and above which they continue to move away from each other. For in-plane interaction the theory allows 
prediction of the outcome of a collision through the inwards escape angle, i.e. whether the solitons 
fuse or cross. We find an analytical expression determining the inwards escape angle using Gaussian 
approximations for the solitons. The theory is verified numerically. 



Stable self-guided laser beams or optical bright spa- 
tial solitons are of substantial interest in basic physics 1 1 ] 
and for technical applications, such as inducing fixed 2] 
and dynamically reconfigurable waveguides Several 
types of spatial solitons have been demonstrated experi- 
mentally, including ID Kerr solitons and ID and 2D 
(number of transverse dimensions) solitons in saturable 
||, photorefractive ||, and x media [7j. Even in- 
coherent solitons, excitable by a light bulb, have been 
demonstrated in photorefractive media ||]. All solitons 
exist when diffraction is balanced by the nonlinear self- 
focusing effect. In bulk Kerr media the self-focusing ef- 
fect dominates and leads to collapse of both coherent and 
incoherent 2D solitons || , their existence requiring an ef- 
fectively saturable nonlinearity. 

An intriguing feature of solitons is their particle like 
behavior during collision. In ID Kerr media collisions 
are fully elastic due to integrability of the ID nonlin- 
ear Schrodinger (NLS) equation flCfl . In contrast, sat- 
urable, photorefractive, and x^ media are described by 
non-integrable equations and soliton collisions are there- 
fore inelastic, displaying both fusion (Fig. [l]:A), cross- 
ing (Fig. |l]:B), repulsion, and annihilation, additional 
solitons can be generated in a fission-type process Jll) , 
and solitons can even spiral around each other JO]. All 
processes depend strongly on the relative phase and have 
been demonstrated experimentally (see [Q and references 
therein). 

Complex waveguide structures can be generated by 
soliton interaction, such as directional couplers |T^| , but 
their efficient implementation requires a detailed under- 
standing of the nature of soliton collisions. Snyder and 
Sheppard predicted the outcome of collisions of ID soli- 
tons in saturable media by comparing the collision angle 
with the critical angle for total internal reflection in an 
equivalent waveguide JLlj . Except for this work most the- 
ories are based on the variational approach, which require 
the solitons to be far apart and breaks down at collision. 

Here we focus on x^ materials [ [j"5| , which are more 
general than the simpler cubic Kerr and saturable media 



in the sense that dependent on the phase-mismatch be- 
tween the fundamental and second-harmonic (SH) waves 
the nonlinearity can be both purely quadratic (close to 
phase-matching) and effectively cubic (for a large phase- 
mismatch). Spatial solitons in x ^ materials do not mod- 
ify the refractive index, and consist of one (type I) or two 
(type II) fundamental fields resonantly coupled to a SH. 
The x materials are of significant interest to photon- 
ics due to the strong and fast nonlinearities they can 
provide through cascading fll6| . Furthermore, soliton in- 
duced waveguides in photorefractives can have a strong 
X^- 2 ' nonlinearity, which can be used for second-harmonic 
generation (SHG) In). 
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FIG. 1. A,B: Planar collision between two Gaussians, 
P = 48.3 and j3 — 0. C,D: Outwards launched exact soli- 
tons, P — 122.4 and j3 = 5. The launch angles are a x — 58° 
(A), a x = 62° (B), a x = 5.4° (C), and a x = 5.7° (D). 

Fusion and crossing (Fig. |:A,B) of spatial x (2) soli- 
tons has been demonstrated numerically p8 and ex- 
perimentally (19), and fission of ID type I solitons was 
demonstrated analytically and numerically in the large 
phase-mismatch limit approximately described by the 
NLS equation pflj. However, the x^ system is more gen- 
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eral and complex than the saturable NLS equation and so 
far variational theories were only able to predict critical 
launch angles and relative phases separating regimes of 
collision and no collision fl21|| . Elegant non-planar effec- 
tive particle theories predicted the absence of spiraling 
type I solitons, but still required weakly overlapping soli- 
tons p^| . In this letter we extend the effective particle 
approach to arbitrary launch angles and present the first 
theory able to correctly predict the outcome of collisions 
between 2D type I solitons in x^ media. 

We consider beam propagation under type-I SHG con- 
ditions in lossless bulk x^ materials. Neglecting walk- 
off the system of normalized dynamical equations for 
the slowly varying envelope of the fundamental wave, 
E x = Bi(f), and its SH, E 2 = E 2 (f), are f|| 



id z E 1 + IVl-Bi + E\E 2 = 0, (la) 
id z E 2 + \V\E 2 - (3E 2 + E\ = 0. (lb) 



Here f=(x,y, z), z is the propagation variable, and 
Vj_ = d 2 + dy accounts for diffraction in the transverse 
rj_ = (x, y)-plane. The normalized phase mismatch is 
(3 — ld(2ki — where Id is the diffraction length of the 
fundamental and k± and k 2 are the wave numbers of the 
fundamental and SH, respectively. The system (|l|) can 
be derived from the Lagrangian density 



C = 21m {E x d z ED + Im (E 2 d z E* 2 



P\E 2 \ 2 



IV^I 2 + ±|V ± £ 2 | 2 - Re (E*E 2 ) 



(2) 



and conserves the power P— J ({E^ 2 +\E 2 \ 2 )dr± and mo- 
mentum M= J Im{£jVi£i + ^E 2 W±E 2 }df±, where we 
have defined J df± = J J_ dxdy. 



The system (]1|) is known [p4|]25| to have a one- 
parameter family of radially symmetric bright 2D soli- 
tons of the form Ei(r) — V(r; A) exp(iAz) and E 2 (r) = 
W(r; X) exp(i2Xz) where A > max(0; — (3/2) is the in- 
ternal soliton parameter and r = \J x 2 + y 2 . We have 
found this family numerically, using a standard relax- 
ation method, and approximately, using the variational 
approach |25|] with Gaussians profiles (V, W) — (V g , W g ), 



tan(a a , ) j / ) are the initial transverse velocities correspond- 
ing to the launch angles a XjV with respect to the z-axis. 

We substitute a field composed of two weakly over- 
lapping solitons (V^- l \W^) into the Lagrangian density 
(^). We then follow the procedure outlined in j2^| and 
allow the solitons to vary adiabatically through a slow 
variation of the soliton parameters with Z = ez being 
the slow propagation variable. To first order in e C 1 
the result is a set of dynamical equations governing the 
collective coordinates x^ l \ yW, and 4>^ l \ being the center 
positions along the x and y-axis and accumulated phase 
of soliton i=\,2, respectively. We can express the new 
coordinates as x®(z) = f* is x l) (Z')dZ' + x { j\ y^(z) = 
f* 4> (Z')dZ' + „W , and M {z) = J* A W (z')dZ' + <#> , 
where subscript denotes initial values. 

At this point one traditionally simplifies the system by 
assuming the velocities to be small (d z x^ ~ e, d z y^ ~ 
e), i.e. the solitons propagate almost in parallel. How- 
ever, we are interested in velocities that can be consid- 
erable, so instead we assume symmetric interaction be- 
tween in-phase solitons with initially identical profiles, 
A = AW and P = PW j and equal but opposite velocities, 
~ Vx \ — —Vx\- Without loss of generality we set 
c (1) - 



-xf ] > and y = y§> 



-yi 2) > 0. Sym- 



X = Xq 

metry is conserved and the two sets of collective coordi- 
nates degenerate to one, X = x^ 1 ' = — x^ 2 \Y = yW = 
— y' 2 ' . In cylindrical coordinates with R = V 'X 2 + Y 2 we 
can then reduce the dynamical equations to the Euler- 
Lagrange equation of the effective Lagrangian 
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L(R,R) = -PR 2 -U eS (R,R), 



for the single coordinate R. The effective potential 



(5) 



(6) 



Vg = aicxp(-r 2 /fr), W g = a 2 exp(-r 2 /6). 



(3) 



is composed of the classical centrifugal barrier, where 
Co = (XY — YX) 2 = {xQVy + yov x ) 2 is constant be- 
cause of conservation of angular momentum, and of the 
interaction integral 



Here ai = a 2 [2(Xb - a 2 = |(A + b^ 1 ), and 

b = [1 + (12A 2 + 8A/3 + /3 2 )i/(2X + /3)]/2A. Because 
the system is Galilean invariant we can apply a gauge 
transformation to find moving solitons. Thus the general 
three parameter soliton family (V, W) is given by 

V(x - v x z, y - v y z\ A, v x , v y ) = 

V(r; X — \v 2 x - \v 2 y) exp[-i(v x x + v y y)], (4a) 
W{x - v x z, y - v y z; A, v x , v y ) = 

j".: ^Vy) exp[-2i(v x x + v v y)}, (4b) 



where (V, W) are either the exact soliton profiles (V s , W s ) 
found numerically or (Vg,W g ) given by (||). i/ x>y = 



u = W 



V {1) W (2) cos(20) + 2W {1) V {2) cos(0) 



df± 
(7) 



where <p = 2v x x + 2v y y. We note that, strictly speaking, 
U is only a quasi-classical potential since it depends on 
the velocities (in contrast to the potential used in p2[). 

We have now established a picture of an effective par- 
ticle moving in a potential, U e g, with the kinetic energy 

^"kin = \PR 2 — 0- F° r small velocities the potential (]?]) 
has the shape of a well and hence represents an attrac- 
tive force. In the general case of non-planar interaction, 
Co 7^ 0, the centrifugal barrier is always repulsive and 
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goes to infinity at R — 0. This does not necessarily rule 
out fusion since also the velocities go to infinity because 
of conservation of Co . The centrifugal barrier also creates 
a local minimum in the effective potential (still assuming 
small velocities) which suggests that spiraling configura- 
tions may exist. In general, however, we cannot expect 
our model to yield correct physical results in the vicin- 
ity of R — since it violates the assumption of weakly 
overlapping solitons. In fact fusion has been observed 
numerically, but stable spiraling configurations have not 
been found p^ |. 

Here we shall not discuss the qualitatively different 
regimes. Rather we are interested in quantitative pre- 
dictions of escape velocities. For solitons launched with 
outwards velocities we will all ways be able to theoret- 
ically predict the escape velocity. On the other hand, 
a consistent theory for the determination of the inwards 
escape velocity only exists for in-plane interaction, when 
the classical centrifugal barrier vanishes, i.e. Co = 0. 

We first determine the outwards escape angle. For 
simplicity we focus on in-plane interaction with yo = 
v y = Co = wherefore R = \X\ and R = v x . In this 
case the effective particle either escapes the potential, 
E tot = ^kin + ^pot > 0, or is trapped by it, E lot < 0, 
and the escape velocity, i/ c , is given by the relation 



v 2 c P = U(x ,v e ). 



(8) 



Unfortunately we are not able to express the interaction 
integral U(xq, v x ) in terms of analytical functions and we 
cannot use Gaussians, since the Gaussian tale asymptotic 
is different from that of the exact soliton. It is however 
trivial to solve (||) numerically and in Fig. || we have plot- 
ted the outwards escape angle, given by a c = Arctan(V c ), 
versus the phase mismatch. The initial beam width and 
separation are kept constant to ensure a weak overlap of 
the soliton tails at all phase mismatches. 




FIG. 2. Outwards escape angle in degrees (solid curve) for 
FWHM—1 and xq— 1.5. Numerical experiments where exact 
solitons fused (•) and where they escaped (o). The dashed 
line shows the initial soliton power, P. 

The simulations were performed with numerically 
found exact solitons as initial conditions. They confirm 
the accuracy of the escape angle predicted by Eq. (||). 
We found the minimum of about 3° around /3 = — 7 to 
be global. Note that the angles are expected to be small, 



since the initial overlap and hence the attractive force be- 
tween the solitons is weak. As an example of the dynam- 
ics we show in Fig. |l|:C,D the outcome of the experiments 
with /3=5 from Fig. |[ Only the fundamental waves are 
shown, the evolution of the SH waves being qualitatively 
the same. 

Now considering solitons launched towards each other 
we first elaborate on the effective particle picture. If we 
assume the solitons to be initially far apart this corre- 
sponds to the effective particle experiencing essentially 
no potential. Even a small launch angle should then re- 
sult in a positive total energy and enable the effective 
particle to cross the bottom of the potential and escape 
towards infinity, corresponding to soliton crossing. This 
is off course not the correct physical picture, since our 
system is not integrable and thus in reality the collision 
is not elastic. There is transfer of energy into internal 
soliton modes and shedding of energy as radiation. 

In a different picture we assume that the soliton profiles 
do not change before the point of collision. This seems 
reasonable when comparing the characteristic length of 
slow adiabatic change with the relatively short interac- 
tion distance occurring for considerable velocities. In this 
case we can treat the interaction as if the solitons were 
launched on top of each other (xo=0) corresponding to 
the effective particle being launched in the bottom of 
the potential, where it experiences the maximum bar- 
rier. Then the relation determining the escape velocity 
becomes v^P — 217(0,^) rather than (||). The interac- 
tion integral no longer depends on the asymptotic tales 
but on the entire profiles and hence we can apply the 
Gaussian approximation (Q). The general transcendental 
equation for the inwards escape angle is then given by 



2 A6+1 
b2Xb- 1 



4 !, 2 



which for f3=0 simplifies to 

13 = : v c = 0.23 x VP, 



(9) 



(10) 



in terms of the power. In the large phase- mismatch cas- 
cading limit, iJ3 A) where the nonlinearity is effectively 
cubic, Eq. (0) simplifies to 



/3> A 



W5 



(11) 



We remark that this approach is equivalent to finding the 
critical angle of total internal reflection for a waveguide 
[ [Til l. However, since beam propagation in quadratic me- 
dia does not induce changes in the refractive index the 
method used in JlJ] is not applicable to this case. 

In Fig. ^ we have summarized the results for exact 
phase matching, [3=0, and plotted the predicted inwards 
escape angle, a c — Arctan(f c ), versus soliton power, 
both for v c given by Jl0| ) and for v c found with exact 
solitons. The curves are close and the simple square 
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root dependency on the power excellently predict the es- 
cape angle. In the experiments we used Gaussians as 
initial conditions. These were launched with a distance 
of 2:eo=10 between them, ensuring practically zero initial 
overlap. In Fig. [l]:A,B we show examples of experiments 
with (3 = and a power of P = 48.3. We note that for 
the exact soliton initial conditions we observed even bet- 
ter agreement than with Gaussians. Close to the escape 
angle all the power is shed as radiation and thus JlO| ) 
serves as an accurate prediction of soliton annihilation 
(Fig.|:B). 



inwards escape angle represents the first analytical pre- 
diction of the outcome of a soliton collision. We have ver- 
ified the analytical expressions numerically using Gaus- 
sian approximations and observed excellent agreement. 

We acknowledge support from the Danish Technical 
Research Council under Talent Grant No. 56-00-0355. 
Much of the numerical work was carried out at Centre 
de Computacio i Comunicacions de Catalunya 
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FIG. 3. Inwards escape angle in degrees versus soliton 
power calculated analytically with Gaussians (solid) and with 
exact solitons (dashed) for P~0 and xo—5. Numerical exper- 
iments where Gaussians crossed (o) and fused (•). 

We also investigated the cases of non-zero mismatches, 
focusing on (3 = ±3. In these regimes there is a power 
threshold for soliton excitation and the collisions are of a 
much more complex nature than for perfect phase match- 
ing, where solitons exist at all powers. For relatively low 
powers not far above the threshold the collisions mostly 
resulted in destruction of the solitons (Fig. ||:A). The ex- 
planation of this phenomenon is that too much power is 
shed as radiation in the collision and hence the resulting 
beams diffract because they do not carry sufficient power 
to form solitons. For higher powers the predicted escape 
angles were reasonably close to the observations. 




FIG. 4. A: Crossing followed by diffraction (/3 = —3 and 
P = 45.3). B: Annihilation (/3 = and P = 48.3). 



In conclusion we have developed a theoretical descrip- 
tion that should hold for systems with all types of local 



nonlinearitics. In particular we have studied bulk \ 



(2) 

media and determined analytical expressions for the es- 
cape angles when the centrifugal barrier vanishes. This 
happens in the two in-plane cases of outwards and in- 
wards launched solitons. The simple expression for the 
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